The Critical State as a Steady-State solution of Granular Solid Hydrodynamics 
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The approach to the critical state - the transition from partially elastic to perfectly plastic behav- 
ior - is considered the most characteristic of granular phenomena in soil mechanics. By identifying 
the critical state as the steady-state solution of the elastic strain, and presenting the main results 
as transparent, analytic expressions, the physics of this important phenomenon is clarified. 

PACS numbers: 81.40.Lm, 46.05.-|-b, 83.60.La 



Applying a constant shear rate to an elastic body, the 
shear stress will monotonically increase - until the point 
of breakage. Granular media are different, as they can 
maintain a uniform state, continually deforming yet of 
constant stress - which has, for given density, a unique 
value independent of the shear rate. This perfectly plas- 
tic state is hailed as a hallmark of granular behavior in 
soil mechanics [Hll], taught in introductory courses, and 
called critical, to convey the idea that it is on the verge 
of stability. All engineering theories are constructed to 
yield a realistic rendering of this state and the approach 
to it - though the physics behind the phenomenon re- 
mains murky. For instance, why is the critical state at 
all stable, and the approach to it smooth? Or, why is it 
rate-independent, but uniform dense flow is not? 

In what follows, we shall first give a brief introduc- 
tion to the present thinking on the critical state; then 
present granular solid hydrodynamics (gsh), a continual 
mechanical theory that, starting from a few simple ideas 
about the basic physics of sand, deduce many aspects of 
its behavior - including the critical state that is identi- 
fied as a steady-state solution of GSH. The main ideas of 
GSH are: • Grains with enduring contacts are elastically 
deformed, which gives rise to static shear stresses. • The 
deformation is slowly lost when grains rattle and jiggle, 
because they briefly lose contact with one another. As a 
consequence, a constant shear rate not only increases the 
deformation, as in any elastic medium, but also decreases 
it, because grains will jiggle when sheared past one an- 
other. A steady state exists in which both processes bal- 
ance, such that the deformation remains constant over 
time - as does the stress. This is the critical state and 
the reason for its perfect plasticity. 

The present thinking is best understood by considering 
a very simple elasto-plastic model, which patches linear 
elasticity with perfect plasticity: Plotting shear stress as 
versus shear strain Ss, the first portion, given by as ^ Ss, 
is elastic; the second, as — ac — constant, is perfectly 
plastic. As the critical stress CTc is equal to the maximal 
elastic one, the critical state is taken to be on the verge 
of stability. Consequently, the point at which elasticity 
switches to plasticity is labeled the yield point. 



More realistically, one needs to include the density de- 
pendence. Assuming cylindrical symmetry, the stress has 
two independent elements, ai — a2 and a^, in the sys- 
tem in which it is diagonal. The approach to the criti- 
cal state is then more varied. The experiments are usu- 
ally performed keeping one stress element constant, say 
ai, [or equivalently P = ^{ai + a2 + a^)], while apply- 
ing a slowly increasing strain, say £3. Employing a tri- 
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FIG. 1: Seemingly a textbook illustration of the critical state, 
this is the result of a calculation employing GSH, specifically 
the solution of Eqs pi7p . plotting the stress q = 0-3 — cti and 
the void ratio e against the strain £3 in tria:xial tests (cylinder 
axis along 3), at given a\ and strain rate zzji, for an initially 
dense and loose sample, respectively. Insets are variations of 
the critical stress qc and density pc with a\. See the text for 
details, and the caption of Fig [2] for the parameter values. 
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axial apparatus, one measures the second independent 
stress element and the density p. With q = (J^ — <Ji 
and the void ratio e = Pg/ p — 1 (where Pg is the den- 
sity of the grains), the results are frequently presented 
as 5(63), e(£3). Provided no shear zones are formed, the 
findings are as rendered in Fig. [TJ {!) q increases mono- 
tonically in loose systems, but displays a maximum in 
dense ones. Both approach asymptotically the station- 
ary value qc- (2) Dense systems dilate, loose systems 
contract, until a universal Cc (or pc) is reached. (3) All 
stress-strain curves are rate-independent, retaining their 
form for whatever strain rate. (4) The friction angle tpc, 
defined as -s/Stant/Sc = qc/P, is essentially constant, in- 
dependent of a I or P, implying qc ^ oi ~ P. The steady 
state, characterized by qc and Cc for given cri or P, is re- 
ferred to as the critical state, see for more details, 
and [J] for a presentation catered to physicists. 

There are many engineering theories capable of a re- 
alistic account of these observations. A mathematically 
elegant and comparatively simple one is the hypoplastic 
theory 0, 0] that starts from a postulated relation be- 

(1) 



tween the stress Uij and strain rate 



dij — Hijkl Vki — MjVs, 

where Vg = ^/vifUiJ > 0. (The compressional rate, 
vu ~ 0, small in the present context, is neglected; the dot 
in dij implies an appropriate objective time derivative.) 
The tensors Hijki,Aij are taken as functions of (yij,p, 
specified using experimental data. The critical state is 
given by <jij — 0. Nevertheless, being descriptive rather 
than explaining, the hypoplastic model - same as vari- 
ous elasto-plastic ones - is not a complete theory: As the 
important features, especially rate-independence and the 
so-called incremental nonlinearity, are put in via Eq ([IJ , 
there is no way to understand them. The hypoplastic 
model is also narrow: None of the rate-dependent phe- 
nomena, say dense fiow 5], or the transition to elastic 
behavior are accounted for. (The latter is needed for 
describing sound propagation [6| and static stress distri- 
butions.) More fundamentally, hypoplasticity does not 
specify a complete set of state variables, and lacks ener- 
getic and entropic considerations. 

In contrast, starting from a few ideas about the physics 
of sand and derived from general principles, including 
energy conservation and thermodynamics, GSH is a full- 
fledged theory that aims to model granular behavior in 
all its facets [3] . It was first developed to calculate static 
stress distribution for various geometries, including sand 
piles, silos, and point load, achieving results in agreement 
with observation [8| . GSH was then generalized to dy- 
namic situations [9|, producing response envelopes strik- 
ingly similar to those from hypoplasticity. At the same 
time, it clarified when rate-independence holds, and why 
incremental nonlinearity applies. Two recent papers are 
on elastic waves_llO| , showing the quantitative agreement 
to experiment t6l|, and on GSH presented from the point 



of view of soil mechanics ll| . Preprints on compaction, 
dense flow, fluidization and jamming have been submit- 
ted, that on shear band is being prepared. In this Letter, 
we consider the critical state. 

GSH consists of flve conservation laws: for the energy 
w, mass p, and momentum pvi, an evolution equation for 
the elastic strain Uij, and balance equations for two en- 
tropies, or equivalently, two temperatures. Both the true 
and the granular temperature, T and Tg, are necessary, 
because granular media sustain a two-stage irreversibil- 
ity: Macroscopic energy, kinetic and elastic, dissipates 
into mesoscopic, inter-granular degrees of freedom, which 
in turn degrades into microscopic, inner-granular degrees 
of freedom. The first are granular jiggling, quantified by 
Tg] the second mostly phonons, quantified by T. The 
elastic strain Uij is defined as the portion of the strain 
Eij that deforms the grains and leads to reversible stor- 
age of elastic energy. The plastic rest, 



accounts 



for rolling and sliding. Because the energy depends on 
Uij alone, we take it as a state variable, while excluding 
the total strain and the plastic one. This is a crucial step 
that enables us to retain many useful features of elastic- 
ity, especially an explicit expression for the stress. 

From the equations of GSH as derived in Q, we need 
three: for Tg, Uij, and the Cauchy stress aij{p,Uki)- 
Accounting for the first-stage irreversibility, the balance 
equation for Tg has a similar structure as that for the 
true temperature, with the production of granular en- 
tropy given as Rg — rjgVg—jTg. The first term, preceded 
by the viscosity rjg, is positive. It describes how the shear 
rate Vij jiggles grains, converting macroscopic kinetic en- 
ergy into Tg. The second term, negative, accounts for in- 
elastic collisions, as a result of which Tg diminishes with 
the rate 7. In steady-state, quickly arrived at (say after 
lO"'^ s), we have Rg = 0, or 



\fhh- 



(2) 



Dividing Uij into A = —uu and u^j = Uij — ^^u uSg, the 
evolution equation for Uij takes the form 0, 

5tA + (1 - a)vu ~ a2u°jVij = -3A/ti = -SAiT^A, 

- (1 - a)v,, = -4 /r = - AT<,4 , (3) 

If Tg is finite, grains jiggle and briefly lose contact with 
one another, during which their unstressed form will be 
partially restored. Macroscopically, this shows up as 



a slow relaxation of A and 



with relaxation rates 



that grow with Tg, and vanish for Tg = 0, we there- 
fore take l/r = XTg and 1/ri = XiTg. Eliminating 
Tg using Eq we have XTg = Avs, iXiTg = AiVs, 
with A = A -^Z 773/7, Ai = 3Ai 1/773/7. Because the relax- 
ation of Uij should come to a halt at the random closed- 
packed density pep, where the system turns elastic, we 
take A,Ai ~ {pcp — p). a,a2 > are Onsager coeffi- 
cients. The first reduces the portion of Vij that deforms 
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the grains and changes the elastic strain u^ ; the second 
quantifies dilatancy, the cross effect typical of granular 
media that a shear flow v^^ leads to a compressional rate, 
(9f A. For the reason discussed below Eq ([3]), we also take 
Ck2 ^ {pcp — p)- The Cauchy stress is 



,{uke,p) = (1 - a)TTij - a2U°.Tru/3, 



(4) 



where TTij{uki,p) = ~dw/duij, w being the elastic en- 
ergy. Viscous terms (large only for the high shear rates 
typical of dense flows) are neglected. That the same a, (X2 
appear here is a result of the Onsager reciprocity relation. 
The term ~ may usually be neglected, as it is an order 
higher in u^, which is always smaller than 10"^. 

Eqs p|4p with Tg eliminated are algebraically simple 
and physically transparent. It is therefore instructive 
that they lead directly to Eq ([T]), the postulated hy- 
poplastic relation, and provide expressions for Hijki , Aij , 
see . Although these tensors appear different from the 
ones used in modern hypoplastic models, the calculated 
response envelopes (ie. closed stress or strain curves) are 
strikingly similar, see Q. This agreement indicates the 
appropriateness of our starting points, especially the ex- 
clusion of the plastic strain as a state variable. 

Next, we evaluate the steady-state solution for the elas- 
tic strain, given by taking dtA, dtu^j = in Eqs As- 
suming a stationary shear flow, = VijVij, vu — 0, and 
denoting u 



^ ^ u°,u°^, we find 

''J '■J 



1 - a 



I - a Ac _ a2_ FY ^ j-^-j 
A ' Uc 3Ai y % Ai ' 



and u^ij/ug — Vij/vg. Reusing the parameter values of Q: 
A = 10^, Ai = 30 for p/p^p = 0.96, a = 0.8 (which are 
the high-Tg limit of a and ijg/j), we obtain Uc = 2 x lO"'^, 
Ac/uc — 0.1 q;2. As the friction angle ipc is a function 
of Ac/uc alone, see Eq (|9]), it determines the value of 
a2 ■ The approach to the steady state is given by solving 
Eqs ^ for Us(t), A(i), at constant p,Vs = Sg/t, with the 
initial conditions: A = Aq, Us = 0. The solution is 



Us{t) - Uc 

A{t) ~ A, 



-Ae, 



a2Uc e 



(6) 



-Ae, 



A- Ai 



Ao 



AAc 

A- Ai 



-Aie, 



As we shall soon see, these simple expressions, giving 
Wc, Ac, Ws(t), A(t) as functions of p and £s = Vgt, are 
a complete account of the the critical state and the 
approach to it. They describe exponential decays of 
Us{t), A{t) to Uc, Ac, and are shear rate- independent, be- 
cause they depend on Eg, not Vg- Given p,A,Ug, the stress 
Uij is known via Eq (|4]). 

The behavior rendered in Fig [T] appears more com- 
plicated, which stems from two facts of more technical 
nature: First, P or ai is usually held constant in stead 
of the density. As A,Us change with time, the density 
compensates to maintain P(p, A, Ug). And along with p. 



the quantities a2,A,Ai, all functions of p, also change 
with time. [In addition, there is a change of time scale 
from the term vu ~ —dtp/ p in the first of Eqs 

Second, the stress is measured, not the elastic strain. 
To calculate <Jij{t) employing Eq we need an ex- 
pression for the elastic energy w. The one we have con- 
sistently employed [1, Q is: w = SVA (fA^ + w^/^) 
where B, ^ are two elastic coefficients. We fix ^ = 
5/3, independent of the density, and take B{p)/Bo = 
[{P - P*)/{Pcp - p)]°'^^, with 9p* = Pep + 2Q{pip - Pep), 
and Bo around 8 GPa for river sand, 7 GPa for glass 
beads. The associated pressure P = k<^ee, and shear 



stress ag = (J^ja^. (with cr° J (Jg = vP^/ug) are 



P = 



(l-a)6Ai-5 [1 + ^2/(2^^2)] 

<jg = {l~a)2ug^B/i, 
Pjog = {^/2)A/ug + {l/A)ug/A, 

where the critical pressure and shear stress are 



Pc = P(p, Ac 



(Tc = ag{p, Ac,Wc). 



(7) 
(8) 
(9) 

(10) 



The expressions Eqs (I5l6l7l8ll0p . evaluated for constant 
(Ti or P, are remarkably similar to textbook illustra- 
tions of the critical state, see Fig [T] and [2] But all 
the well-taught features of these curves are now easy to 
understand: We first note that while Us{t) always in- 
creases monotonically, A{t) decreases monotonically only 
for f2 = [Ao- AcA/(A - Ai)] > 0, as /i = a2Wc/(A-Ai) 
is always positive. We also note that for given P, a posi- 
tive /2 means stronger initial compression, corresponding 
to a lower density or higher void ratio, while /2 < sig- 
nifies weaker compression and lower void ratio. At the 
beginning, the faster relaxation of /i dominates, so A, e 
always decrease, irrespective of the void ratio e. After 
/i has run its course, A,e go on decreasing if /2 > 0, 
but switch to increasing if /2 < 0, displaying respectively 
the so-called contractancy and dilatancy, until the critical 
state, A = Ac, e — Cc, is reached. 

The shear stress ag always increases first with Ug, until 
Ug — Uc is reached. The subsequent behavior depends on 
what A does [cf. Eq Q and the discussion on Coulomb 
yield below, consider only Ug/A < -\/2f]: ag keeps grow- 
ing if A decreases [loose case, /2 > 0], but decreases, 
displaying a peak, if A grows [dense case, /2 < 0]. 

The friction angle tennpc = Qc/V^P — ag/\/2P is ob- 
served to be essentially independent of the density, or 
the pressure [H [3]. This is also accounted for by the 
above results, because the ratio CTs/P depends only on 
Ac/uc — Q;2/Ai, see Eqs ([5l and we have taken 
0^2, Ai ^ Pep — p. There are two types of density depen- 
dence here: the sensitive one via pcp — p, and the weaker 
one via p. Neglecting the latter, a2/Ai is a constant. 
Circumstances are similar if P is the control parameter, 
because P ^ B ^ {pcp — p)""'^^ for constant Uij, so 
changing P only changes pcp — p. 
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The relation of to Coulomb yield is instructive. 
Given Eq ([9]), the ratio P/as has the minimal value 
y/ij2, at u^/A = V2?- It is larger for m^/A < ^/2^, 
and unstable for Ug/A > (as the energy w is then 
concave, see [8|). This may be identified with (the Druck- 
Prager version of the) Coulomb yield - an energetic in- 
stability that ensures that no elastic solution exists for 
P/ag < or Ms/A > \/2^. Remarkably, this insta- 

bility is determined by the elasticity coefficient ^, while 
the friction angle is given by transport coefficients, espe- 
cially a2 ■ There is no a priori link between both, as they 
are based on different physics. So the critical state may 
be stable and observable, or not. For the first case, we 
have Mc/Ac = Ai/a2 < implying tan ipc < 

or (fc < 38° for ^ = 5/3. (Taking as before Ai = 30, this 
is equivalent to a2 > 3^/30.) Clearly, if critical is meant 
to imply marginal stability, it is a misnomer - energetic 
stability is essential for the existence of the critical state. 

The yield point discussed in the introduction is also a 
dubious concept. Lacking an obvious choice in the upper 
two curves of Figdl one frequently takes it as given by q's 
maximum for dense sand - the vague rational being the 
fact that in an elastic medium, the positivity of the stiff- 
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FIG. 2: Variation of (a) the shear stress as, (b) density p, 
(c) granular temperature Tg, (d) volumeric deformation rate 
Vkk, (e) bulk strain A, (f) deviatoric elastic strain 113—111, with 
the total shear strain Ss, of a initially dense (po = 0.97pcp) and 
loose (po = 0.95pcp) sample, for fixed P and axial deformation 
rate. Initially, both as and Tg vanish. Parameters: Bo — 
7GPa, ^ = 5/3, pep = 0.85pcp, Q2 = 750(1 - p/pcp), a = 0.7, 
A = 2850(1 - p/pcp), A/Ai = 3.3, and in SI units: 7 = 6 x 
lO^Tg, rig = 0.15 Tg. (Note A = 114, aa = 30 for p = 0.96pcp.) 
The steep linear curve (red) is the purely elastic case. 



ness coefficient, dq/de^lai = da^/de^ > 0, is required by 
energetic stability, so a negative slope implies instability 
and is never observed. However, since Uij is the state 
variable and not Sij, stability only requires da^/du^ to 
be positive, freeing da^/de^ to be negative, as observed. 



In spite of the energetic stability of the critical state, 
the yield surface may of course still be breached at some 
point. Shear bands (considered in a forthcoming paper) 
will then appear, destroying the uniformity of the system. 
This is most likely to happen around g's maximum for 
the densest sand, because Uc meets the smallest A there. 



Summary: Sheared granular media, if uniform, will 
approach the continually deforming critical state that is 
widely believed to be on the margin of stability. Employ- 
ing GSH, a broadly applicable granular theory, we find the 
critical state well accounted for, with all facets transpar- 
ently explained, and given by the steady-state solution 
of the elastic strain that is both continuous and stable. 
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